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Abstract: A cavity QED implementation of the non-adiabatic holonomic 
quantum computation in decoherence-free subspaces is proposed with 
nitrogen-vacancy centers coupled commonly to the whispering-gallery 
mode of a microsphere cavity, where a universal set of quantum gates 
can be realized on the qubits. In our implementation, with the assistant 
of the appropriate driving fields, the quantum evolution is insensitive to 
the cavity field state, which is only virtually excited. The implemented 
non-adiabatic holonomies, utilizing optical transitions in the A type of 
three-level configuration of the nitrogen-vacancy centers, can be used to 
construct a universal set of quantum gates on the encoded logical qubits. 
Therefore, our scheme opens up the possibility of realizing universal holo¬ 
nomic quantum computation with cavity assisted interaction on solid-state 
spins characterized by long coherence times. 
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1. Introduction 


Quantum computation is an attractive field on precision control of quantum systems. But quan¬ 
tum systems are inevitably influenced by the decoherence effect induced from the environment, 
which stands in the way of physical implementation of quantum computation. In order to over¬ 
come this difficulty, there are many strategies to correct or avoid errors during the implemen¬ 
tation of the gate operations. The geometric phases null and their non-Abelian extensions, 
quantum holonomies EH), accompanying in evolutions of quantum systems, unveil important 
geometric structures in the description of the dynamics of quantum states. Geometric phases 
depend only on global geometric properties of the evolution paths so that they are largely in¬ 
sensitive to many types of local noise. Based on this distinct merit, holonomic quantum com¬ 
putation (HQC), first proposed by Zanardi and Rasetti 0, has emerged as a promising strategy 
to implement universal quantum computation in a robust way 00. 

It is well known that the geometric phases, either Abelian or non-Abelian, consist of both 
adiabatic and non-adiabatic parts. In adiabatic evolution, it requires that quantum evolutions 
fulhll the famous adiabatic condition, and thus quantum gates constructed in this way generally 
have very slow speeds. This is unacceptable as the time needed for an adiabatic quantum gate 
maybe on the order of the coherence time of the used quantum two-level systems in typical 
systems Colin]. Therefore, one needs to consult quantum gates based on non-adiabatic evolu¬ 
tion cniini, which does not require the adiabatic condition so that the gate speeds will only 
depend on the merits of the employed quantum systems. In this case, the built-in fault-tolerance 
of the non-adiabatic geometric phases provide a more practical way in implementing quantum 
computation. Therefore, many renewed efforts have recently been given in this direction both 
theoretically and experimentally ICMHI- 

Meanwhile, as an another promising way to avoid the effect of decoherence, decoherence- 
free subspace (DFS) can suppress the collective dephasing noise caused by the interaction be¬ 
tween quantum systems and their environment ||T5l420l . Therefore, many efforts have been paid 
to combine the HQC with DFS encoding ll2TI425ll . In this way, one can consolidate the best of 
the two quantum computation strategies, i.e., resilient against environment induced collective 
decoherence of the DFS approach and the operational robustness of the HQC against local 
noises. 

Here, we propose a non-adiabatic HQC (NHQC) scheme with DFS encoding based on 
nitrogen-vacancy (NV) centers coupled commonly to the whispering-gallery mode (WGM) of 
a fused-silica microsphere cavity. The NV center system is considered as a promising candidate 
for physical implementation of quantum computation, due to its sufficiently long electronic spin 
lifetime as well as the possibility of coherent manipulation even at room temperature 02611271 . 
Since the electronic spins of the NV centers can be well initialized and manipulated in an 
optical way, quantum gates acting on the single-spin state can be obtained with very high ef- 
hciently 02811290 . Based on the symmetry structure of the system-environment interaction, we 
propose an encoding method of the logical qubits in the DFS, where any logical qubit state that 
undergoes a cyclic evolution will end up with a new state in the same subspace, without going 
out of the subspace. Moreover, based on the numerical simulation under realistic conditions, 
we show that our NHQC scheme in DFS can realize a universal set of quantum gates with high 
hdelity. 

2. Setup and effective Hamiltonian 

The setup we consider and the energy level configuration of the NV centers are schematically 
shown in Fig. [T] In the cavity, the lowest-order WGM, corresponding to the light traveling 
around the equator of the microsphere, offers exceptional mode properties for reaching strong 
light-matter coupling. As shown in the Fig. N NV centers, from separate diamond nano- 
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Fig. 1. (a) Schematic setup of the fused-silica microsphere cavity, where N identical NV 
centers (red dots) in diamond nanocrystals are equidistantly attached around the equator of 
the cavity, (b) Level diagram for a NV center, with g (Q.l) is the coupling strength between 
NV center and the WGM (laser pulse). We encode a physical qubit in the subspace spanned 
by the states of = 0 and = — 1. 


crystals, are strongly coupled to the WGM of a microsphere cavity ll^lSTI . The NV centers 
we considered can be modeled as three-level systems, as shown in Fig. 1(b), where the states 
\^A,ms = O) and \^A,ms = —l) are encoded as our qubit states |0) and |1), respectively. The 
state \^E^ms = O) is labeled by \e) and leave aside the metastable *A state, which has not yet 
been fully understood ll^ . In our implementation, the optical transition |0) \e) and |1) —^ \e) 

(with transition frequencies (Ogo and Cdgi) are coupled by the WGM with frequency (Og and a 
classical laser filed with frequency col 13^ . respectively. Both coupling is far-off resonant from 
their transition frequencies so that the \e) state can be adiabatically eliminated. Note that there 
is a degenerate state \^A,ms = l) for the qubit state \^A,ms = —l), these two degenerated states 
can be selectively addressed by the polarization of the optical field lf3^ . The NV centers are 
fixed and the distance of two NV centers is much larger than the wavelength of the WGM, so 
that each driven laser field, with frequency cOl and initial phase (pj, can interact individually 
with an NV center and the direct coupling among the NV centers can be negligible. Then the 
interaction of the whole quantum system, in units of h = 1, can be written as 


Ho 


Hs 


0)ca'^a + Y^C0i\i){i\, 


N r 


f=i 


( 1 ) 


where (a) is the creation (annihilation) operator of the WGM of the cavity, O); is the frequency 
of /th energy level of the identical NV centers with i G {0,1,e}, t]/ {^lj) is the coupling 
strength between jth NV center and the cavity (laser). In the interaction picture with respect to 
the free Hamiltonian Hq, under the rotating-wave approximation (RWA), the interaction can be 
written as 

Hi=Y Sj + H.c.) , (2) 

)=i 

where we have assumed rjj = G for simplicity, the effective cavity assisted interaction strength 
are gj = GD.lj ( { ) with A = cUgo - ®c, (7^^ = 11) ,• (01, ar = 10)(11, and 5±j = ±{(Oc- 
coio-colj). 

When two laser fields are applied to a pair of NV centers (for example the mth and the nth), 












one can obtain the effective Hamiltonian, under the condition of 5± 3> gmn, as 
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where (p^„ = (p,„ — (p„. Neglecting the level shift terms, which can be compensated by using 
additional lasers II341I35L the effective Hamiltonian between the two NV centers reduces to 


17 _ 6mn / i<p„, 

0± 


C^in CTy, 


H.c.) 


(4) 


With this Hamiltonian, we next show how a universal set of non-adiabatic holonomic gates can 
be implemented. 


3. Non-adiabatic holonomic One-qubit logical gates 

Now we turn to construct universal single logical qubit operations for the NHQC in DFS. It is 
noted that a logical qubit consists of two physical qubits is not sufficient for a dephasing envi¬ 
ronment We here utilize three physical qubits to encode a logical qubit. The interaction 
between three physical qubits (NV centers) and the dephasing environment can be described by 
the interaction Hamiltonian, H[ = (^E, where <7? is the collective dephasing op¬ 

erator and E is an arbitrary environment operator. In this case, there exist a three-dimensional 
DFS of 

5i={|100), 1001), 1010)}, 

where the computational basis, i.e., the logical qubit states, are encoded as |0 )l = 1100), 11)/, = 
|001), and |ai) = |010) is used as a third ancillary state. 

In order to implement the single logical qubit operation, we apply operation Unit) = 
exp[—/ Jq* Hndt], which can be obtained from Eq. (01) for an operate time of Ti, to the NV 
centers 1 and 2 by tuning laser beams with the same detuning 5+ from their respectively transi¬ 
tion and phase difference (p = (pi — (pi- Meanwhile, we also apply t/ 23 ('Ci) = exp[—i H 2 idt] 
on NV centers 2 and 3 by tuning laser beams with the same detuning 5 and identical phase. 
Then, in the DFS the effective Hamiltonian reads 

Hi = Ai ^sin^e'‘'’|ai)z,(0| -cos^|ai)L(l|^ +H.c., (5) 

where the effective Rabi frequency Ai = \/|gi 2 |^ +1^23|^/5 with 5 = |5±| and the phase 
0 = 2arctan (|gi 2 p/|^ 23 p) can be tuned by the amplitude of the incident lasers. This Hamilto¬ 
nian, consisting of a WGM and three physical qubits (three NV centers), is in the A type with 
ancillary logical state \a\)L at the top while the logical qubit states |0)l and | 1)l at the bottom. 
In the dressed state representation, the two degenerate states of Eq. dSl) are 

\d)L = cos^|0)L-f sin^e'‘'’|l)z„ 

I^)z. = sin^e^"P|0)L-cos^|l)L. (6) 

It is obviously that the dark state \d)L decouples from the ’bright’ state \b)L and the excited 

state |ai)L, while the ’bright’ state \b)L couples to the excited state |ai)i with the effective Rabi 
frequency Ai. 

As a result, the evolution operator U\ (Ti ) = exp(—i H\dt) can realize the holonomic gates 
with following two conditions. The first one is that controlling the operation time to meet Ai Ti = 
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Fig. 2. Qubit states population and fidelity under single-qubit operator U\{6,(p) with initial 
state |0)/,, the red (dot) line represents the fidelity, (a) 9 = Ttjl and (p = 0. (b) 6 = njA and 
(p = 0. 


n. This condition will ensure the states evolving in the subspace \'^fi{t))L = Ui{t)\i)L{i = d,b) 
undergo the cyclic evolution |tp, (Ti))L = | V/, (0))l. Furthermore, the relation | v/y (f)) = 

0 {i,j = d,b) is always hold, which means the evolution is pure geometric and the parallel- 
transport condition is natural satisfied. Therefore, we can obtain the non-adiabatic holonomic 
single qubit gates in the space spanned by {|0)l, |1)l} as 


Ui{0,(p) 


cos 9 sin0e \ 
sinOe'^ —cos 9 J’ 


(7) 


where 9 and <p can be chosen by tuning the amplitude and relative phase of laser beams, and 
thus a set of universal single qubit gates can be realized non-adiabatically. For examples, we 
can implement a Hadamard gate with t/i(;r/4,0) and a phase gate as Ui(7l/2,7l/4)Ui(7r/2,0). 

Inevitably, the implementation process will suffer from decoherence. Considering the main 
decoherence effect, the collective relaxation rate 7 , dephasing rate 7 ^ of NV centers and the 
decay rate K of the cavity, we simulate the performance of our scheme under realistic conditions 
with the Lindblad master equation ll^ 

p =—i[Hp\p] +-[K^{a) + Y^{S ) + ( 8 ) 

where is the Hamiltonian in the form of Eq. (|2]i for the case of = 3 with the compensation 
of the Stark shift, p is the density operator, and J/f {s^Y) = 2£/pj!/^ — sYp — 

pxY^xY is the Lindblad operator. In our simulation, we have used the following conservative 
set of experimental parameters. The NV centers are located near the microcavity surface and 
the maximum coupling between an NV center and the cavity could be G = 27r x 1 GHz with 
the mode volume of Vm = lOOpnP IIJTI . For 5 ^ A, g ~ 2GQ.l/A = 27r x 50 MHz with Q.l = 
271 X 500 MHz, A = 27r X 8 GHz and 5 = 27r x 1 GHz to fulfill the condition of 5 ^ g. The qubit 
relaxation and dephasing rates are estimated to be 7 = 7 ^ w 2;r x 4kHz lIMl . The cavity decay 
rate is K = (odQ ~ 2;r x 0.5MHz with Q = 10^ ll^ . Assume that the logical qubit is initially 
prepared in |0 )l state while the cavity is in the vacuum state, the time-depend state populations 
and fidelity under the X and Hadamard gates are depicted in Fig. 12 a) and|2b) with the fidelity 
to be 99.5 and 99.6, respectively. Note that, in the simulation, we use the Hamiltonian in Eq. 














Fig. 3. Maximum fidelity as the function of 0 in unit of TZ with line (circle) the identical 
(different) decoherence of qubits under gate operations (a) Ui{tz/2, 0), (b) U\ (?t/4,0). 


dU, and thus the obtained high fidelity also verifies the validity of the effective Hamiltonian 
in Eqs. Q and (|4|i. In addition, we assuming the collective dephasing of all the qubits, the 
violation of which will introduce very small infidelity to the gate operations, as shown in Fig. 
|3a) and[3^b) for the two exemplified gates f/i(7r/2,0) and U\{n/A,0) respectively, where we 
have investigated different cases by choosing different 0 in the initial state of 11/) = cos 0 |O)l + 
sin0|l)L. In the simulation, we have choose different decoherence of the physical qubits as 
7i = 7<|>i = 0 - 87 . 72 = 7<|>2 = 7 and 73 = 7^3 = 1 . 27 , respectively. 

4. Non-adiabatic holonomic two-qubit logical gates 

The implementing of the holonomic one-qubit logical gates can be scalable to two-qubit sce¬ 
nario straightforwardly. For two logical qubits interacting collectively with the dephasing envi¬ 
ronment, there will exist a six-dimensional DFS 


S2 = { 1100100 ), 1100001 ), 1001100 ), 1001001 ), 1101000 ), 1000101 )}, 


where the former and latter three physical qubits encode the first and second logical qubits, 
respectively. Then we can encode the logical qubit states same as that of the single qubit case, 
i.e., |00)l= 1100100), |01)l = |100001), |10)l = |001100), and |ll)i = 1001001). Meanwhile, 
1 ^ 2 ) = 1101000) and |a 3 ) = |000101) are used as ancillary states. 

In order to implement the logical two-qubit operations, we apply the interaction in Fq. (|4|l on 
the NV centers 3 and 4, by tuning laser beams with detuning 5+ and phase difference <j) = (p^ — 
(p 4 , for the operation time of T 2 , i.e., 1 / 34 ( 72 ) = exp[—) jQ^H 34 dt]. Meanwhile, we also apply 
— exp[—! H^sdt] on the NV centers 3 and 6 by tuning laser beams with detuning d 
and same phase. The total effective Hamiltonian then reads 


H2 = ^2 


sinye''^(|a2)L(00| 


|a3)L(ll|)-cosy(|a2)L(01| + |fl3)z,(10|) 


-FH.c., 


(9) 


where the effective Rabi frequency X 2 = •\/|g 34 |^ + /8 with 5 = |5±| and the phase 

tJ = 2arctan (|g 34 p/|g 36 p) can be tuned by the amplitude of the incident laser. Obviously, 
the effective Hamiltonian can be decomposed to two commuting parts as H 2 = ^liHa + Hj,) 
















Fig. 4. Qubit states population and fidelity under two-qubit operator (/2(^/4,0), the red 
line represents the fidelity, (a) Initial state of qubit |00)r,. (b) Initial state of qubit lOl)/,. 


with 


Ha = sin |a2)z.(00| - cos y |a2)z.(011 + H.c., 

Hb = sinye''^|fl 3 )i(ll| -cos Y|a3)z,(10| +H.C.. (10) 


Therefore, exp [—/ H 2 dt] = exp[—/7 i://a]exp[—under the K pulse criterion A 2 T 2 = TZ, 
which acts nontrivially on the computational subspace {|00 )l, |01)l} and {|10)L,|ll)i}, re¬ 
spectively. That is, H 2 , consisting of a WGM and six NV centers, effectively reduces to two 
A-like Hamiltonian with Ha {Hb) acts on the top state \a 2 )L ( 1 ^ 3 ) 1 ) and the bottom states |00)i 
and|01)L (|10 )l and |11 )l)- Analogous to the single-qubit gate case, the holonomic two-qubit 
logical gate in the subspace {|00)£, lOl)/., | 10 )l, |1 1)l} can be written as 


U2{^,^) 


/ COS J? 

sin 
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V 0 


sintJe 
— cos 
0 
0 


0 0 \ 

0 0 
— cost? sintJe^''^ 
sin d-e'^ cos t? / 


( 11 ) 


In general, we can realize nontrivial two-qubit holonomic logical gate in DFS by controlling 
the t? and (j) separately, that is, adjusting the amplitude and phase of two laser beams, respec¬ 
tively. For example, acting U 2 {'&,(l)) on logical qubit 1 and 2 and then acting Ui{9,(p) on logical 
qubit 2 with t? = 0 = 7r/4 and (j) = (p = n/2, we can realize a CNOT gate in view of Eqs. (I7]i 
and (fTTll . In addition, choosing the experimentally achievable parameters similar as in the single 
qubit case, numerical simulation of the populations and the fidelity of t/2(?i/4,0) operation are 
shown in Fig.Ua) and|4fb) with different initial state |00 )l and |01 )l, respectively. The fidelity 
can reach about 99.5 and 98.7. 


5. Conclusion 

In summary, we have put forward a universal set of non-adiabatic holonomic gates in DFS by 
using NV centers coupled to the WGM of a cavity. By controlling the amplitude and relative 
phase of the driving lasers, we can realize arbitrary single-qubit and two-qubit operations. Nu¬ 
merical simulation shows that our scheme is stable to deviation of corresponding experimental 























parameters and insensitive to decoherence such as collective noises and local noises. Ultrahigh 
quality factor of the cavity and the exceptional spin properties of the NV centers make our 
scheme a promising candidate in experimental implementation of NHQC in DFS with high 
gate fidelity and short operation time. 
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